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The self-organization of magnetospheric plasma is brought about by inward diffusion of magne¬ 
tized particles. Not only creating a density gradient toward the center of a dipole magnetic field, 
the inward diffusion also accelerates particles and provides a planetary radiation belt with high 
energy particles. Here, we report the first experimental observation of a ‘laboratory radiation belt’ 
created in the Ring Trap 1 (RT-1) device. By spectroscopic measurement, we found an appreciable 
anisotropy in the ion temperature, proving the betatron acceleration mechanism which heats parti¬ 
cles in the perpendicular direction with respect to the magnetic field when particles move inward. 
The energy balance model including the heating mechanism explains the observed ion temperature 
profile. 


I. INTRODUCTION 

Magnetospheres are natural plasma confinement de¬ 
vices ubiquitous in the universe. Yet, their creation 
mechanism remains to be fully understood. There must 
be a spontaneous mechanism that transports particles 
towards the center of a dipole magnetic field [1, 2]. How¬ 
ever, such inward diffusion seemingly violates the en¬ 
tropy principle, because it creates a gradient instead of 
flattening the density distribution. This challenge may be 
overcome by considering adiabatic invariants; these are 
the magnetic moment, the bounce action and the longi¬ 
tudinal angular momentum. Typically, the frequencies 
of the corresponding periodic motions are separated by 
three orders of magnitude. The third invariant is most 
fragile, and its violation (i.e. the change of angular mo¬ 
mentum) gives rise to the radial transport of particles. 
When the diffusion of magnetized particles is constrained 
by the remaining adiabatic invariants, the particle den¬ 
sity tends to distribute uniformly on ‘magnetic coordi¬ 
nates’ rather than on Euclidean coordinates. Therefore, 
the homogeneous density on the former turns out to be 
inhomogeneous on the latter [3-5]. 

The same mechanism plays another interesting role, 
i.e., the acceleration of particles to produce a radiation 
belt [6-10]. Conserving the remaining adiabatic invari¬ 
ants, the kinetic energy of a particle increases as it moves 
inwards. An increase in cyclotron frequency (maintain¬ 
ing the magnetic moment constant) results in betatron 
acceleration of the velocity in the direction perpendic¬ 
ular to the magnetic field, and an increase in bounce 
frequency (maintaining the bounce action constant) re¬ 
sults in Fermi acceleration of the parallel velocity [2] (the 
electrons are further accelerated to an ultra-relativistic 
regime by whistler waves [11, 12]). Since the former is 
stronger than the latter, the inward diffusion results in 
temperature anisotropy [2]. Numerous satellite observa¬ 
tions have provided evidence of anisotropic temperatures 
in planetary radiation belts. The anisotropic electron 
temperatures are consistent with the estimates from in¬ 
ward diffusion heating [13-15]; ions also have anisotropic 


temperatures in the radiation belt of the Earth [16, 17] 
as well as that of Saturn [18]. 

In addition to these theoretical and observational stud¬ 
ies, the inward diffusion was experimentally verified re¬ 
cently. On the Ring Trap 1 (RT-1) device which is a 
‘laboratory magnetosphere’ simulated by a levitated su¬ 
perconducting magnet [19, 20], the peaked electron den¬ 
sity profiles were observed and proved to originate from 
inward diffusion [21]. Similar density profiles were ob¬ 
served in the Levitated Dipole Experiment (LDX), which 
is another dipole confinement system [22]. On RT-1, dis¬ 
tinct proof of inward diffusion was provided by produc¬ 
ing a non-neutral (pure electron) plasma [23]; by action- 
conserving acceleration, electrons diffuse into the central 
region where the electric potential is higher than their 
initial energies (such particles absorb energy from fluctu¬ 
ations, which drive radial motion). Here, we report the 
first experimental observation of the other outcome of the 
inward diffusion, i.e. the particle acceleration. Specif¬ 
ically, we observed the ion temperature anisotropy and 
identified its generation mechanism as betatron acceler¬ 
ation. This finding means that we succeeded to demon¬ 
strate ‘laboratory radiation belt’. 


II. EXPERIMENTAL SET-UP 

The experiment was performed on RT-1 device. Fig¬ 
ure 1 shows the layout of the device, together with asso¬ 
ciated plasma images (soft X-ray image and an electron 
density distribution reconstructed from interferometry). 
Plasma is produced by electron cyclotron resonance heat¬ 
ing (ECH) with an 8.2 GHz microwave (maximum power 
of 50 kW) and a discharge duration of approximately 1 s. 
The plasma contains high-temperature (10 ^ 50 keV) 
electrons as well as low temperature (< 100 eV) elec¬ 
trons; typically, the former occupies about half of the 
total population [24]. Electrons are confined in a radially 
elongated and vertically thin region, resembling plane¬ 
tary radiation belts [21, 25]; see Fig. 1(a). From quasi- 
neutrality, ion density distribution is also considered to 
be peaked on the equatorial plane (z = 0). In the inward 
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diffusion process, a concomitant heating mechanism must 
be in effect to produce high-energy electrons. However, 
electrons are also heated by ECH (cyclotron resonance 
occurs just in the vicinity of the levitated ring magnet). 
It is, therefore, difficult to separate the betatron and 
Fermi accelerations in the total energy balance of elec¬ 
trons. As a viable alternatives, ions may thus be used 
as an appropriate target for analyzing potential heating 
mechanisms. 

Spatial distribution of ion temperature was measured 
by two sets of Doppler spectroscopies; one scanned the 
plasma on a horizontal plane, and the other on a vertical 
plane (see Fig. 1). The lines of sight of the former ranged 
from r = 375 mm to 878 mm, and those of the latter 
ranged from r = 375 mm to 808 mm on the z = 0 axis (r 
and 2 : are the radial and vertical coordinates, respectively, 
r = 375 mm is the surface of the superconducting mag¬ 
net, r = 1000 mm is the vacuum chamber wall, and z = 0 
is the mid plane of the magnetic dipole). Hell transition 
(468.565 nm) was used to measure the Doppler broaden¬ 
ing. The horizontal chords measured the perpendicular 
temperature (T±) with respect to the ambient magnetic 
field, while the vertical chords measured the mixture of 
T± and the parallel temperature (T||). The local tem¬ 
perature profiles of T± and T\\ were reconstructed from 
the line integrated spectrum data, (see Appendix for the 
detailed algorithm). 



FIG. 1: (a) The plan (vertical front view) of the device. The 
superconducting magnet is levitated by the lifting magnet 
placed on the top of the vacuum chamber. The lines of sight 
of the vertical spectroscopic measurement cover the greyed 
vertical plane. On the left-hand side of the figure, the elec¬ 
tron density profile (measured by interferometers [21, 25]) is 
plotted and the superimposed soft X-ray image (the interfer¬ 
ometry measurement and the soft X-ray image are from dif¬ 
ferent shots) is shown, (b) The plan (top view) of the device. 
The lines of sight of the horizontal spectroscopic measurement 
cover the greyed horizontal plane. The reflectometer (3.9 GHz 
microwave) measures the density fluctuation near the cutoff 
density which is ~ 1.9 x 10^^m“^ (light blue). 


III. OBSERVATION OF ION TEMPERATURE 
ANISOTROPY 


Figure 2(a) and 2(b) shows the two dimensional profiles 
of T|| and of the He+ ion in helium plasma. Because 
of the mirror effect, T|| decreases and T±_ increases with 
closer distance to the magnetic poles. The radial tem¬ 
perature profile of the He+ ion on the equatorial plane 
is plotted in Fig. 2(c), showing an evident difference be¬ 
tween Tx and T\\. 
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FIG. 2: Two dimensional profiles of (a) parallel temperature 
(T||) and (b) perpendicular temperature (T^) ofHe^ ions. T±_ 
is plotted in a region where the mirror ratio (the local mag¬ 
netic field normalized by the minimum magnetic field along 
each field line) is less than 2. (c) The radial temperature 

profile on the equatorial plane {z — 0). The translucent re¬ 
gions are the error bands estimated by the covariance of the 
reconstruction fitting. 


Figure 3 shows the dependencies of line averaged ion 
temperatures and anisotropy on ECH power and neu¬ 
tral helium particle density at two different radial points 
(r ^510mm and r ^710mm). The neutral helium par¬ 
ticle density is estimated as = p/ksT — where p 
is the filling gas pressure, ks is the Boltzmann constant, 
T = 300 K is the room temperature and is the electron 
density measured by the interferometers. From Fig. 3(a) 
and 3(b), both the temperatures and anisotropy increase 
as ECH power increases. From Fig. 3(c) and 3(d) both 
the temperatures and anisotropy decreases as the neutral 
helium particle density increases. 
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FIG. 3: Dependence of (a) the line averaged temperature of 
He^ and (b) the anisotropy on the ECH power, and the de¬ 
pendence of (c) the line averaged temperature and (d) the 
anisotropy on the neutral helium gas density. 


IV. ENERGY TRANSPORT MODEL 
INCLUDING BETATRON ACCELERATION 


nisms depends on the driving mechanism of the inward 
diffusion. 

Figured shows the frequency spectrum of the density 
fluctuations detected by a reflectometer. A flat spec¬ 
trum extending up to around 10 kHz was observed. The 
bounce frequency of a He+ ion is typically around 15 kHz 
(assuming Ty ^ 10 eV). On the other hand, the ion cy¬ 
clotron frequency is typically 360 kHz (for = 20 eV 
and B = 150 G); thus, the betatron acceleration is 
deemed to be ideally in effect. Hence, it was concluded 
that in this instance, the Fermi acceleration is at most 
very weak in comparison with the betatron acceleration. 
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In order to distinguish the effect of betatron acceler¬ 
ation as the preferential heating mechanism for T^, the 
energy balance in the He+ ions was examined. By conser¬ 
vation of the magnetic moment, must increase upon 
displacement towards the central region as the magnetic 
field strength increases. From conservation of the first 
adiabatic invariant, it may be estimated that AT^/Tj^ = 
AH/H, where B is the magnetic field strength. Let Ar 
denotes the radial displacement by an infinitesimal time 
interval At. The heating power by betatron acceleration 
is given by 


betatron 


ATj_ 

At 


Tl_Vr 


dlnS 

dr 


( 1 ) 


where Vr = Ar/At is the speed of inward diffusion. 

A close parallel mechanism is the Fermi acceleration 
which heats Ty via an increase in the bounce frequency. 
In a point-dipole magnetic field (which may approximate 
a planetary magnetic field), it is known that betatron 
acceleration yields much stronger heating than Fermi ac¬ 
celeration [2]. In the magnetic field of RT-1, on the other 
hand, both of them may yield almost same powers if the 
first and second actions are equally conserved. But if 
the constancy of the second action is broken upon in¬ 
ward diffusion, only the betatron acceleration will be in 
effect. The relative importance of these heating mecha- 


FIG. 4: Time evolution of the frequency spectrum measured 
by the reflectometer for the total lifetime of the plasma (EGH 
power is injected from 0.5s through 1.5s). The density cut¬ 
off resides at r ~ 640 mm. The inset is the clip between 
t = 1.001584 s and t = 1.063193 s. The fluctuation has a 
flat spectrum up to around 30 kHz (which ranges beyond the 
bounce frequency). However, the fluctuations in the range of 
the ion cyclotron frequency are very weak. 


Other competitive processes are (1) the thermal equili¬ 
bration of He+ ions and the cold component of electrons, 
(2) the thermal equilibration of He+ ions and He+ ions 
immediately after ionization from neutral helium parti¬ 
cles which have isotropic temperature Tq, (3) the loss of 
He+ ions by charge exchange with neutral particles and 
(4) the isotropization of and T\ \. The effect of hot elec¬ 
trons may be neglected because the thermal equilibration 
time is on the order of 10^ s, which is considerably longer 
than the lifetime of the plasma. In this instance, no as¬ 
sumption is made of any ‘anomalous transport’; thus, the 
only energy loss mechanism at play is charge exchange 
loss. 

Taking into account the aforementioned processes, a 
one-dimensional model of energy transport on the radial 
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coordinate r may be considered: 
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where d/dt is the temporal derivative to be evaluated 
on each volume element co-moving with the plasma on 
a Lagrangian coordinate; for a stationary state, we may 
evaluate d/dt = Vrd/dr with an inward diffusion speed 
Vr. All coefficients of equations ( 2 ) and (3) are evaluated 
by experimentally measured local plasma parameters. 

The electron-ion thermal equilibration time (rei) is a 
function of the ion temperature 10 eV), the cold elec¬ 
tron density (ng ^ 10 ^^ m“^ measured by the interferom¬ 
eters) and the cold electron temperature Tg [26]. Figure 5 
shows the spatial profile of Tg measured by line integrated 
He I line ratios (728.1/706.5 and 667.8/728.1 nm) [27] 
for different plasma condition from that of Fig. 2. Un¬ 
like the ion temperature profile, Tg is higher near to the 
dipole field magnet. We only have the single spatial point 
measurement Tg = 30.3 eV at r = 795 mm for the same 
condition as Fig. 2. Therefore, we infer that the cold 
electron temperature is distributed as 30 < Tg < 40 eV 
in the plasma for the present study. Finally, Tgi was de¬ 
termined to be ^ 0.1s. The thermal equilibration time 
among He+ ions (ru) is estimated as ^ 1 ms by the ion 
temperature and density ng). The ion isotropization 
time (riso) is on the order of 1 ms, as estimated by the 
ion temperatures (T^ ^ 20 eV and T|| ^ 10 eV) [26]. 


neutral gas and is the ionization rate coefficient. 

The neutral helium temperature was found to be almost 
spatially homogeneous with 0.7 eV by the Doppler broad¬ 
ening of the He I (471.315 nm) line as shown in Fig. 6 . 
Using rig of the shot in Fig. 2 and 30 < Tg < 40 eV, A 
is evaluated as ^ 1 m near the dipole field magnet and 
^ 10 m at the chamber wall, implying that the neutral 
gas density is spatially flat. In addition, the absolute 
value of the neutral helium density must be evaluated. 
We solved the rate equation of ionization, recombina¬ 
tion and charge exchange processes of helium, varying 
neutral helium gas density. The stationary solutions are 
shown in Fig. 7. Since adequate light emission was ob¬ 
served from both the He I and Hell lines, the neutral 
gas density was determined to be roughly in the range 
of O.lrig ^ iin ^ lOrig. If rin is too small, He^+ dom¬ 
inates and He I line would be scarce. Or if rin is too 
large. Hell line would disappear (since He^+ does not 
emit light, we have yet to determine how much of He^+ 
is present). Therefore, the neutral helium gas density is 
homogeneously rin ^ 10 ^^ m“^. Tgx was determined to 
be ^ 1 ms. 
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FIG. 5: The spatial profile of Te measured by line inte¬ 
grated He I line ratios (728.1/706.5 and 667.8/728.1 nm). The 
plasma condition is different from that of Fig. 2. 

Let us estimate neutral helium gas density profile to 
evaluate the time constant Tgx which scales the time of 
charge transfer between the He+ ion and the neutral he¬ 
lium gas. The mean free path of the neutral particle is 
A = Vu/tIq (cr'^)ion’ where Vn is the thermal velocity of the 


FIG. 6: Radial profile of the line averaged temperature of the 
neutral helium particles evaluated from the Doppler broaden¬ 
ing of the He I line spectrum (471.315 nm) on the horizontal 
line of sight. 


The inward diffusion speed Vr is required to evalu¬ 
ate Tbetatron- This was measured in the initial forma¬ 
tion phase. At the beginning of the discharge, the den¬ 
sity has a broad distribution, but subsequently it begins 
to concentrate into the central region, thereby creating 
a clump of particles. Figure 8 (a) shows the time evo¬ 
lution of the line integrated electron density measured 
by the three chord interferometers around the onset of 
the plasma discharge. Two dimensional density profile 
of rig is reconstructed at the two different time points 
t = 0.98300 s and t = 0.98348 s. The electron density 
evolves in accordance with the continuity equation with 
the inward diffusion velocity V : 


drie 


1 ^ 

r dt 


(rrigW) 


dz 


{rieVz) = 0 . 


( 4 ) 
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FIG. 7: Ionization ratio of helium as a function of neutral 
particle density, obtained by solving the rate equation of ion¬ 
ization, recombination and charge exchange processes. The 
electron temperature was set to 30 eV from the measurement. 


We assume V as 


V = V^o 



a 


ivv^r 


The parameters and a are optimized so that the 
solution of (4) with the initial condition as dX t = 
0.98300 s in Fig. 8 becomes closer to ng at t = 0.98348 s. 
Choosing appropriate parameters, the density evolution 
is replicated well as shown in Fig. 8(a), thus the speed W 
that transports the particles is deduced. In the stationary 
phase, we do not have a direct measurement of Vr^ but 
we assume that the same estimate applies. Figure 8(b) 
shows optimized Vr on the equatorial plane which will be 
used to evaluate Pbetatron- 



FIG. 8: (a) Time evolution of the line integrated elec¬ 

tron density measured by the interferometers at each chords 
around the onset of the plasma discharge. The bullets at 
t = 0.98348 s are the solution of (4) with the parameters 
(Wo, ol) = (3000, —7) and the initial condition as the bul¬ 
lets at t = 0.98300 s. (b) Radial distribution of the estimated 
inward diffusion speed Vr on the equatorial plane. 


V. NUMERICAL SIUMLATION OF THE 
ENERGY TRANSPORT MODEL 

Equations (2) and (3) are both ultimately solved with 
the aforementioned parameters. An initial condition at 
r = 1000 mm (edge of the plasma) is assumed; using 
the measurement, T± = E|| =8 eV. This initial tem¬ 
perature may be equivalent to the temperature of ions 
immediately after ionization from neutral helium parti¬ 
cles. Therefore we set Tq = 8 eV. The solution is shown in 
Fig. 9(a) and it turns out that the model explains the ex¬ 
perimental profiles presented in Fig. 2(c). In the first half 
of the inward transport (1000 > r > 700 mm), the be¬ 
tatron acceleration is stronger than the charge exchange 
loss because of the sufficiently fast W. This results in 
the increase in T±. In the last half (700 > r > 375 mm) 
where W is rather smaller, on the other hand, the charge 
exchange loss overcomes the betatron acceleration. Thus, 


T± decreases. Since riso is smaller than the inward dif¬ 
fusion time scale in this region, the isotropization also 
decreases Tl, and at the same time, increases T|| Here 
the exact value of the neutral gas density (which has 
a homogeneous distribution in the plasma) is tuned to 
1 X 10^^ m“^ to match the peak value of Tj_. 

Figure 9(b)-9(e) illustrate the dependencies of temper¬ 
atures and anisotropy calculated from our model on the 
initial temperature at r = 1000 mm and the neutral he¬ 
lium particle density. The results are consistent with 
the experimental scaling in Fig. 3. The increase in the 
ECH power implies an increase in the initial ion tem¬ 
perature through the relaxation between cold electrons; 
then Pbetatron bccomcs larger, resulting in larger T± and 
anisotropy. The increase in T| | is given by the isotropiza¬ 
tion. The decrease in the ion temperature with the in¬ 
crease in neutral helium gas density is consistent with the 
charge exchange loss term in (2) and (3). Since Pbetatron 
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is independent of the neutral gas density, an increase Acknowledgments 

in the charge exchange loss results in small anisotropy. 
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In conclusion, the anisotropic ion temperature ob¬ 
served in the laboratory radiation belt is experimental 
proof of the betatron acceleration concomitant with in¬ 
ward radial diffusion of particles. 
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FIG. 9: (a) The radial profiles of T± and T|| calculated by 
the energy transport model (2)-(3). The result is compared 
with the experimental observation shown in Fig. 2(c). Depen¬ 
dence of (b) the ion temperature and (c) the anisotropy on 
the initial temperature at r = 1000 mm, and the dependence 
of (d) the ion temperature and (e) the anisotropy on the neu¬ 
tral helium gas density. The temperatures and anisotropy are 
evaluated at r = 700 mm. The results are compared with the 
experimental scaling shown in Fig. 3. 
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Appendix A: Algorithm for the reconstruction of 
the local profile 

To reconstruct the local temperature profile from the 
line averaged measured data, the following procedure was 
used. The local profile on the equatorial plane (z = 0) 
is reconstructed from the spectrum data obtained via 
the horizontal chord. From the axisymmetric property, 
the spectrum on the equatorial plane can be written as 
/eq(A, '0(r, 0)) where A is the wavelength, '0(r, z) is the 
inverse of magnetic fiux and r and z are the radial and 
vertical coordinates, respectively. The innermost fiux 


surface contacting the dipole field magnet is denoted as 
and the outermost fiux contacting the vacuum cham¬ 
ber at r = 1000 mm is denoted as Hereafter, the 
subscript eq means the value on the equatorial plane. 
The line integrated spectrum ^eq(A, y) is thus calculated 
from the local spectrum as 

5 'eq(A, y) = 2 f /eq(A, ^eq{'/x^^^))dx, (Al) 

Jo 

where x denotes the auxiliary coordinate along the line 
of sight. We assume 


/eq(A, V^eqW) = < 


Aeqi'lpeqir)) eXp 


A - Ao - Beqii’eqir)) 


Ceq(V’eq(?’)) 


(■^^0 < 'Ipeqir) < i’l) 


{ipl < Ipeqir)) 


where Aq is the central wavelength. The intensity The line of sight passing the equatorial point (^, 0) is 
Aleq(' 0 eq(^)) ^ud the shift Seq(' 0 eq(^)) ^rc assumcd to be given as 


^eq('0eq(^)) — 


{^eq/^0 - - V^eqM)^" 


r=— — -{z-Zc)+rc. (A3) 



Cisi^O + 0 ^ 2^01 
(22 + (23 


The line integrated spectrum ^(A, y) is calculated from 
the local spectrum /(A, r, z) as 


^eqii’eqir)) = bi + b2lpeq + (>3^/'eq> 

with the perpendicular temperature on the equatorial 
plane assumed to be 


7iLeq(t^eq) 


a A'^’eq/'lpO-iy^il-'lpeq/'lplY^ 


TO.5^4 


1 T erf 





tsYo + hYi \ 

t2+t3 ) 


■ (A 2 ) 


The Doppler broadening Ceq('0eq) is determined by 
7±eq('0eq)- The functious ^eq^Seq and T^eq indicate 
that the spectral emission is zero at both boundaries, 
the temperature is zero at -0 = -00 and ti at -0 = t/^i, 
and the toroidal flow speed may be finite at both bound¬ 
aries. The observed spectrum data set for the horizon¬ 
tal chord is written as (A^, yj^ ^eqij)- The parameters 
(ai, a 2 , as, 6 i, 62 , 63 , ti, ^ 2 , ts, t^) are optimized so as to 
minimize |^eqij - ^eq(Ai, %)| for each A^ and yj. 

Next we build a vertical reconstruction algorithm. Us¬ 
ing ^eq and Teq obtained by the horizontal reconstruc¬ 
tion, a two-dimensional profile is reconstructed. The po¬ 
sition of the collimator is denoted as (rc, z^) (Fig. 1(a)). 


In the same manner as the horizontal reconstruction, the 
local spectrum is assumed to be 


/(A, r, z) 


I A{r, z) exp 


fX-XoV 

\C{r, z)) 


{Yo<Y < Yi) 


0 


(■^^1 < Y) 


In the vertical chord, the Doppler shift can be ignored 
because of the absence of the mean poloidal flow in the 
dipole magnetic held configuration. We use the mag¬ 
netic held coordinate to treat the vertical direction, and 
assume A{r, z) to be the power function of the magnetic 
held strength B. 


4(r, B{r, Z)) = AeqiYeq) 


B{r, z) 


B, 


eq 


(A5) 


where B^^ is the magnetic held strength of the same held 
line on the equatorial plane. In the z 7 ^ 0 region, the 
effective temperature (Tgff) causing the Doppler broad¬ 
ening is a combination of and T\\. Defining the angle 



0 between the line of sight and the line normal to the 
magnetic field, we may write 

(Teff)2 = (Tx cosdf + (T|| sin 0)2. (A6) 

In the same manner as (A5), we assume the perpendicu¬ 
lar temperature profile as 

= (A7) 

While T^eq is already determined by the horizontal re¬ 
construction, we need the profile of T|| on the equatorial 
plane. We assume T||eq in a similar manner as T±eq- 


r||eq(V’eq(r)) 


(Ji — o.sq) 


ii’eq/i’O - 1)*"(1 - V>eqM)*" 
{'Ip^/'lpO - 

+0.5(4 [l+e.A 


0.2{'tpi - V’o)/. 



h + h ) 


.(A8) 


In the same manner as (A5) and (A7), we assume 


(^) horizontal • measured — optimized 

(r>“rtical • measured — optimized 



y [mm] 


FIG. 10: Line averaged temperature profile to estimate the 
fitting precision. The bullets are the line averaged tempera¬ 
tures of He^ estimated by the measured spectrum data. The 
solid lines are the line averaged temperatures estimated by 
the optimized g. The translucent regions are the error bands 
determined by multiplying a constant by the covariance of the 
reconstruction fitting. 


7||{r,4B(r, z)) = T||eq(^T_^ (A9) 

We write the observed spectrum data set for the 
vertical chord as (A^, yj^ gij). The parameters 
(a. Cl, C 2 , ti, ^ 2 , ts, ^ 4 ) are optimized to minimize \gij — 
g{Xi, yj)\ for each and yj. A two-dimensional pro¬ 
file of A{r, z), T±(r, z) and T'||(r, z) is then ultimately 
obtained. 

Figure 10 shows the fitting precision of Fig. 2(c). The 
solid lines are the line averaged temperatures estimated 
by geq{X^ y) and ^(A, ^), respectively, with optimized pa¬ 
rameters, while the bullets are the line averaged temper¬ 
atures estimated by the measured data ^eq ij and gij , re¬ 
spectively. The translucent regions are the error bands 
determined by multiplying a constant by the covariance 
of the reconstruction fitting. The corresponding error 
bands for the local profile are illustrated in Fig. 2(c). 
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